Interaction networks within disease-associated GαS variants characterized by an integrative biophysical approach

Activation of G proteins through nucleotide exchange initiates intracellular signaling cascades essential for life processes. Under normal conditions, nucleotide exchange is regulated by the formation of G protein–G protein–coupled receptor complexes. Single point mutations in the Gα subunit of G proteins bypass this interaction, leading to loss of function or constitutive gain of function, which is closely linked with the onset of multiple diseases. Despite the recognized significance of Gα mutations in disease pathology, structural information for most variants is lacking, potentially due to inherent protein dynamics that pose challenges for crystallography. To address this, we leveraged an integrative spectroscopic and computational approach to structurally characterize seven of the most frequently observed and clinically relevant mutations in the stimulatory Gα subunit, GαS. A previously proposed allosteric model of Gα activation linked structural changes in the nucleotide-binding pocket with functionally important changes in interactions between switch regions. We investigated this allosteric connection in GαS by integrating data from variable temperature CD spectroscopy, which measured changes in global protein structure and stability, and molecular dynamics simulations, which observed changes in interaction networks between GαS switch regions. Additionally, saturation-transfer difference NMR spectroscopy was applied to observe changes in nucleotide interactions with residues within the nucleotide binding site. These data have enabled testing of predictions regarding how mutations in GαS result in loss or gain of function and evaluation of proposed structural mechanisms. The integration of experimental and computational data allowed us to propose a more nuanced classification of mechanisms underlying GαS gain-of-function and loss-of-function mutations.

Activation of G proteins through nucleotide exchange initiates intracellular signaling cascades essential for life processes.Under normal conditions, nucleotide exchange is regulated by the formation of G protein-G protein-coupled receptor complexes.Single point mutations in the Ga subunit of G proteins bypass this interaction, leading to loss of function or constitutive gain of function, which is closely linked with the onset of multiple diseases.Despite the recognized significance of Ga mutations in disease pathology, structural information for most variants is lacking, potentially due to inherent protein dynamics that pose challenges for crystallography.To address this, we leveraged an integrative spectroscopic and computational approach to structurally characterize seven of the most frequently observed and clinically relevant mutations in the stimulatory Ga subunit, Ga S .A previously proposed allosteric model of Ga activation linked structural changes in the nucleotide-binding pocket with functionally important changes in interactions between switch regions.We investigated this allosteric connection in Ga S by integrating data from variable temperature CD spectroscopy, which measured changes in global protein structure and stability, and molecular dynamics simulations, which observed changes in interaction networks between Ga S switch regions.Additionally, saturation-transfer difference NMR spectroscopy was applied to observe changes in nucleotide interactions with residues within the nucleotide binding site.These data have enabled testing of predictions regarding how mutations in Ga S result in loss or gain of function and evaluation of proposed structural mechanisms.The integration of experimental and computational data allowed us to propose a more nuanced classification of mechanisms underlying Ga S gain-of-function and loss-of-function mutations.
Signaling at the surface of most eukaryotic cells is triggered by the interaction of agonist-stimulated G protein-coupled receptors (GPCRs) with intracellular G proteins.GPCR-G protein complex formation accelerates the exchange of bound GDP for GTP, dissociating the trimeric G protein complex into a and bg subunits.Both subunits then activate intracellular pathways until GTP is hydrolyzed, enabling the dissociated subunits to reassemble and restart the process (1,2).In typical physiological conditions, nucleotide exchange and the cyclical dissociation and reassociation of G proteins is carefully controlled through interactions of G proteins with partner GPCRs and regulators of G protein signaling proteins.Acquired or inherited single point mutations, or residue-specific modifications, within the Ga subunit of the trimeric G protein subvert this process.Certain single point mutations within the Ga subunit accelerate GDP release or GTP uptake, causing constitutive activation and uncontrolled generation of secondary messengers (3,4).Other mutations prevent GDP release or dissociation of the Ga and Gbg subunits, irreversibly halting secondary messenger production (5).In both cases, circumventing regulation of G protein signaling is closely associated with the onset of a wide range of diseases (6)(7)(8), including inflammatory diseases (9), cardiovascular diseases (10,11), metabolic disorders (12,13), and cancers (14)(15)(16)(17).
The prevalence of disease-associated mutations in multiple Ga subtypes has prompted investigation of the structural mechanisms by which these mutations cause abnormal signaling.Literature data from integrative biophysical and biochemical experiments of the inhibitory Ga protein, Ga i , indicated the presence of a function-related allosteric network connecting the nucleotide-binding pocket and Ga i "switch" regions (18,19), structural features conserved among Ga proteins that undergo conformational changes upon nucleotide exchange (20).The presence of this allosteric network, and the correlation between changes in function with changes in hydrogen bonding between switch regions, also underlies a rational basis for explaining the impacts of disease-associated mutations in Ga S (21).A crystal structure of the Ga S R201C variant bound to GDP revealed a pattern of hydrogen bonding interactions among crucial switch regions similar to those in GTP-bound Ga S (21).This indicated the structure of the GDP-bound R201C variant represented an active conformation, providing a plausible rationale for the gain of function observed for this variant (21).Additionally, the same study proposed that the loss of hydrogen bonding for the GTPbound Ga S variant R228C could explain the loss of function for this disease-associated variant (21).These findings raised the question of the extent to which changes in protein stability and changes in hydrogen bonding within the switch regions correlated with functional changes in other Ga S variants.However, data on these aspects are lacking for most Ga S variants.No structures are available for other Ga S variants, and even for variants that have been structurally characterized, no structural information is available for both GDP and GTP complexes.It thus seemed pertinent to systematically measure these properties across a wider range of Ga S variants.
Motivated by the lack of data measuring stability and hydrogen bonding for most Ga S variants, we leveraged an integrative experimental and computational biophysical approach to characterize the structure and stability of seven of the most prevalent mutations in Ga S and their influence on Ga S structure-function relationships.Our studies were carried out in aqueous solutions, circumventing the need to crystallize the variants, which can be challenging.Using variabletemperature CD spectroscopy, we measured the global structural fold and stability for all variants.Changes in networks of hydrogen bonding and nonhydrogen bonding interactions were observed through molecular dynamics (MD) simulations, which enabled visualization of changes in interactions among the crucial switch regions.To observe the allosteric effects of mutations on the nucleotide-binding pocket, we employed saturation-transfer difference NMR (STD-NMR) spectroscopy in aqueous solutions to examine how mutations impacted the interaction of nucleotides with residues in the nucleotidebinding pocket and provide a structural and dynamic interpretation of this data.These data were compared with predictions made in earlier studies of Ga S , which enabled us to evaluate and revise previously proposed mechanisms of how disease-associated mutations alter Ga S function.Given the occurrence of analogous mutations in other Ga subtypes, the mechanisms explored in this work may extend beyond Ga S .

Prediction of properties of human Ga S variants
Our study compared the biophysical and structural properties of Ga S and 7 Ga S variants containing single amino acid replacements: R201C, G226A, Q227L, R228C, R258A, R265H, and A366S (Fig. 1).These variants were selected because of their disease relevance (Table S1) and because alterations to their adenylate cyclase activity had been previously measured (Table S1).This allowed us to propose a framework for relating the functional impacts of Ga S mutations with the biophysical and structural properties of the variants measured in our experimental and computational data.Selected sites of mutation spanned four important structural regions.R201 is located in switch I; G226, Q227, and R228 are located in switch II; and R258 and R265 are located in switch III (Fig. 1).A366 is located between the b5 and a5 regions adjacent to the nucleotide-binding pocket, approximately 3.5 to 4.0 Å from the imidazole group of the guanine ring system.
Although the stability of Ga S and the variants in the present study had not been previously directly measured via assays that monitor protein unfolding, a framework for predicting their relative stabilities was developed based on earlier structure-function studies.A crystal structure of GDP-bound Ga S [R201C] showed that this mutation altered an allosteric link between the nucleotide-binding pocket and the 3 Ga S switch regions.The structure of GDP-bound Ga S [R201C] revealed an increased number of hydrogen bonds among the switch regions, which made GDP-bound Ga S [R201C] more closely resembled the structure of GTP-bound Ga S (21).An increased number of hydrogen bonds were predicted to correlate with relatively higher stability for GTP-bound Ga S and GDP-bound R201C, therefore correlating an increase in adenylyl cyclase activity for the R201C variant with an increase in the stability of GDP-bound R201C.The correlation between stability, hydrogen bonding networks, and aberrant function suggested a mechanism for the loss of adenylate cyclase activity for other Ga S variants, in particular R228C, would involve an inverse effect.That is, destabilization of the GTPbound conformation was predicted to correlate with loss-ofactivity mutations (21).
Based on the above reasoning, we sought to investigate the extent to which changes in hydrogen bonding networks within the switch regions correlated with global protein stability and to also investigate the extent to which gain-of-function mutations in Ga S shared the property of increased stability for their GDP-bound conformation.Concurrently, we also aimed to determine the extent to which a decrease of hydrogen bonding within the switch regions correlated with decreased global protein stability and to also determine the extent to which loss-of-function mutations shared the property of decreased stability for their GTP-bound conformation (Table 1).We predicted that comparing the stability and hydrogen bonding networks across a greater range of Ga S variants would allow us to understand which variants shared similar mechanisms of gain of function with R201C or loss of function as R228C or whether alternative mechanisms were needed to explain the functional impact of other mutations.
Three of the seven studied mutations are documented to be gain-of-function mutations that increase cAMP accumulation: Ga S [R201C] (21), Ga S [Q227L] (22,23), and Ga S [A366S] (4) (Table S1).For these, we sought to determine to what extent the mechanism proposed for the gain-of-function of Ga S [R201C] could be extended to the additional gain-offunction mutations.Four different loss-of-function mutations have been documented to result in decreased cAMP production: Ga S [G226A] (5), Ga S [R228C] (21), Ga S [R258A] (21) and Ga S [R265H] (21).For these variants, we sought to determine to what extent the mechanism proposed for Ga S [R228C] could be extended to explain loss of function for additional variants.For each of these variants, we next tested these predictions against experimental data measuring their stability (see below).

Preparation of human Ga S and Ga S variants
Human Ga S and all Ga S variants were expressed in Escherichia coli and purified to >95% homogeneity using protocols adapted from earlier studies (21,24); see Experimental procedures and Fig. S1, A and B for additional details.CD spectra of purified Ga S and all Ga S variants showed minima at 210 and 220 nm (Figs.2A and S2), typical of folded proteins containing alpha helical secondary structure, confirming all samples were folded.For Ga S , structural rearrangements in the three switch regions associated with nucleotide exchange result in an increase in intrinsic tryptophan fluorescence (20).To verify protein function, we monitored changes in intrinsic tryptophan fluorescence upon exchange of GDP for the nonhydrolysable GTP analog GTPgS (Fig. S1).Rates of tryptophan fluorescence enhancement for Ga S and Ga S variants were consistent with values reported in prior studies (5,21).

Functional changes caused by mutations partially correlate with predicted changes in protein stability
Previous investigations documented that Ga i1 complexes with GTP or GTP analogs exhibited a $10 C higher thermal denaturation temperature compared with GDP-bound Ga i1 (19,25,26).These observations were interpreted to indicate that Ga i1 adopted different conformations between the GDP and GTP-bound states, with the GTP-bound state exhibiting an increased number of hydrogen bonding among the switch regions giving rise to increased global structural stability (19,25,26).
To measure the stability of the 7 Ga S variants and compare them with the stability of Ga S , we determined their thermal unfolding temperatures using variable temperature CD spectroscopy in the presence of either GDP or GTPgS.To provide an internal control, the study also included "mini-Ga S ," an engineered G protein with a truncation of the alpha helical domain and multiple mutations introduced to facilitate its crystallization in complex with GPCRs (24,27).Mini-Ga S recapitulates certain essential aspects of Ga S binding to GPCRs, however it does not bind nucleotides with the same affinity as Ga S and does not undergo global conformational rearrangements upon nucleotide exchange (24).The unfolding temperature (T m ) of GTPgS-bound Ga S was observed to be $9 C higher than GDP-bound Ga S (Fig. 2, B and C), consistent with earlier findings of Ga i1 (19,25,26).
Among the three variants with gain-of-increase mutations, we observed an increased T m for GDP-bound Ga S [R201C], consistent with the earlier finding of an increased hydrogen bonding network that would result in increased stability (21).However, for both Ga S [Q227L] and Ga S [A366S] we observed a T m value for their GDP-bound conformations to be similar to Ga S (Fig. 2C), in contrast to an initial prediction of a higher Properties and predicted impacts of the studied mutations, including their location, previously reported impact on adenylyl cyclase activity, and the potential impacts on stability and hydrogen bonding network that were tested in this study.See also Table S1 for related properties of these variants.
Disease-associated mutations of G proteins T m .Interestingly, for Ga S [Q227L], we observed similar T m values for both the GDP-bound and GTP-bound samples as compared with Ga S .For Ga S [A366S], we observed a lower T m value for the GTP-bound protein, suggesting a global conformation for the GTP-bound protein that more closely resembles GDP-bound Ga S .Among the four variants with loss-of-function mutations, we observed decreased T m values for Ga S [R228C], consistent with the earlier prediction of decreased stability for the GTPbound conformation of this variant (21).We also observed a decreased T m value for Ga S [G226A], suggesting this variant may share a similar mechanism of loss of function (see Discussion).However, for both Ga S [R258A] and Ga S [R265H], we observed an increased T m value for the GDP-bound conformation and a T m value for their GTP-bound conformation similar to that of Ga S .Applying the same logic used to understand the functional and biophysical impact of the R201C mutation, the higher T m value for GDP-bound Ga S [R258A] and Ga S [R265H] indicates an increased hydrogen bonding network for their GDP-bound conformations, in contrast to the prediction of a correlation between a loss of hydrogen bonding and decreased T m value.
For comparison, we also recorded thermal unfolding data for mini-Ga S .As expected, we observed at most minimal differences in temperature between mini-Ga S in the presence of GDP or GTPgS (Fig. 2), in line with the absence of structural differences observed in the presence of GDP or GTPgS and consistent with reports that claimed mini-Ga S did not respond to the presence of GTP (24).
For the 3 Ga S variants with increased T m values of their GDP-bound conformations, Ga S [R201C], Ga S [R258A], and Ga S [R265H], we sought to understand what role the interaction with bound GDP played in providing increased stability by comparing the T m values of these variants in the presence and absence of GDP.Ga S , Ga S [R201C], Ga S [R258A], and Ga S [R265H] were expressed and isolated without nucleotide added during the purification process.Purification of additional Ga S variants without nucleotide was not attempted or GTPgS (orange curve).B, thermal unfolding of GaS bound to GDP or GTPgS monitored by variable temperature single wavelength CD spectroscopy.Same color scheme as A. C, histogram of Tm values determined from fitting the variable temperature CD data of GaS, mini-GaS, and GaS variants obtained in complex with GDP or GTPgS.Variant names are organized according to those with predicted increased stability for the complex with GDP (cyan) or predicted decreased stability for the complex with GTP (red), relative to GaS.Dashed horizontal lines are located at the mean Tm values for GaS in complexes with GDP or GTPgS.Error bars were determined from the SD of triplicate measurements.
because these experiments focused on comparing Ga S and variants that exhibited increased T m values for their GDPbound states.Apo Ga S exhibited a T m of $62 C, higher than the T m of GDP-bound Ga S and lower than the T m of GTP-bound Ga S (Fig. S3).This suggests apo Ga S adopts a global conformation that differs from both GDP-bound and GTP-bound Ga S .Interestingly, although Ga S [R201C], Ga S [R258A], and Ga S [R265H] all exhibited higher T m values for complexes with GDP, we observed strikingly different behavior for apo preparations of these variants.For apo Ga S [R201C], we observed a lower T m value of around 59 C. Notably, both apo Ga S [R258A] and Ga S [R265H] exhibited higher T m values similar to the T m values in the presence of either nucleotide (Fig. S3).This indicated that the global conformations of Ga S [R258A] and Ga S [R265H] are similar both in the presence and absence of nucleotides but that Ga S [R201C] appears to adopt different conformations in the absence or presence of nucleotides.The implications for these observations on mechanisms of gain of function and loss of function are expanded on in the Discussion.

MD simulations reveal changes in hydrogen bonding that correlate with thermal melting profiles of different Ga S variants
Our framework for understanding the impact of the studied mutations posits that increased stability of Ga S driven by mutations correlated with an increased hydrogen bond network.As discussed earlier, literature data of Ga i1 stability and structure showed a higher T m value for GTP-bound Ga I1 , which was attributed to increased interactions between switch I and switch II and between switch I and switch III upon complex formation of Ga i1 with GTP or GTP analogs (18,19).The structure of GDP-bound Ga S [R201C] revealed an increased hydrogen bond network similar to that of GTPbound Ga S (21).Based on these earlier studies, we hypothesized that Ga S variants with increased T m values for complexes with GDP would exhibit an increased hydrogen bond network interactions, while Ga S variants with decreased T m values for complexes with GTP would exhibit a decreased hydrogen bond network interactions.
To investigate changes in the hydrogen bonding network, we conducted MD simulations on Ga S and several Ga S variants in complexes with GDP and GTP.We focused our efforts on simulating Ga S, the gain-of-function mutation Ga S [Q227L], which showed unexpected thermal melting behavior similar to native Ga S , and the loss-of-function variants Ga S [R228C] and Ga S [R258A], which showed different thermal melting profiles from each other despite that both mutations cause loss of function.
Based on our framework relating protein stability and hydrogen bonding interactions, we predicted that for the lossof-function mutation Ga S [R228C] we should observe a decrease in the interactions within the hydrogen bonding network as compared with Ga S .Analysis of hydrogen bond interactions captured from MD simulations revealed that for Ga S [R228C] we observed significant loss of hydrogen bonds between switch I and switch II and between switch II and switch III for the GTP-bound state as compared with native Ga S (Fig. 3).This is qualitatively consistent with the measured lower T m value for GTP-bound Ga S [R228C] (Fig. 2).We observed the loss of an interaction between Ser205 in switch I and Glu230 in switch II for GTP-bound Ga S [R228C], and the loss of two additional interactions between switch II and switch III that were present for GTP-bound Ga S, between Arg228 and Glu259 and between Arg228 and Glu268 (Fig. 3).Comparing all GTP bound proteins, we observed the fewest number of interactions between switch regions for Ga S [R228C].
In MD simulations of GTP-bound Ga S [R258A], we observed the hydrogen-bonding network to be more similar to that of native Ga S than Ga S [R228C] (Fig. 3).As compared with Ga S [R228C], for Ga S [R258A] we observed several additional interactions between switch I and switch II and between switch II and switch III not observed for Ga S [R228C] (Fig. 3).While not identical to the hydrogen bond network observed for Ga S , this observation is qualitatively consistent with the similar T m values measured for GTP-bound Ga S [R258A] and that of Ga S .For Ga S [R258A] bound to GDP, we also observed a notable increased hydrogen bond network, including increased hydrogen bond interactions between switch II and switch III and the presence of two interactions absent for Ga S , between Gln227 and Glu259 and between Arg228 and Glu259 (Fig. 3), which were absent for GDP-bound Ga S .This increased hydrogen bond network for GDP-bound Ga S [R258A] also appears qualitatively consistent with the increased T m measured for GDP-bound Ga S [R258A] shown in Figure 2. Though both Ga S [R228C] and Ga S [R258A] are loss-offunction mutations, the distinct differences in their hydrogen bonding networks for both GDP-bound and GTP-bound conformations is consistent with their strikingly different stabilities and suggests two different mechanisms underlying their loss of activity (see Discussion).
MD simulations of Ga S [Q227L] showed patterns of hydrogen bonding interactions that were most similar to those of native Ga S among the three variants for both GDP-bound and GTP-bound Ga S [Q227L] (Fig. 3).The only differences observed between Ga S [Q227L] and Ga S were the loss of one interaction between switch II and switch III involving Q227 and E259 for the GTP-bound protein and the addition of one interaction between switch II and switch III for the GDPbound protein.The similarities in hydrogen bond networks between Ga S and Ga S [Q227L] correlate with the differences in measured T m values for Ga S [Q227L] bound to GDP and bound to GTPgS and the T m values for Ga S bound to GDP and bound to GTPgS (Fig. 2).
To complement the analysis of hydrogen bond networks, we also compared the dynamics of each switch region, described by the root mean square fluctuation (RMSF), among Ga S and the Ga S variants considered in this analysis.RMSF values were calculated from averaging over all MD trajectories for each variant and nucleotide combination.Fluctuations were determined as the average backbone deviation relative to the average position of the full trajectory per residue.We observed Disease-associated mutations of G proteins the largest RMSF for Ga S [R228C] for both the GDP and GTPbound forms of this variant (Fig. S4), particularly around residues 225 to 230 of switch II.We also observed increased RMSF values across five residues of switch III for Ga S [R228C] (Fig. S4).This is in line with the decreased hydrogen bonding observed for residues between switch II and switch III of Ga S [R228A] (Fig. 3) and consistent with the observed lower T m value for GDP-bound Ga S [R228C] shown in Figure 2.
Following methods for studying nonbonded interactions among switch regions in Ga i1 (19), the hydrogen bond analysis was further complemented with quantification of nonbonded interaction energies averaged over the MD trajectories, calculated for groups of residues defining the three switch regions (Fig. S5).The total interaction energy between groups of residues in switch I and switch II and between switch II and switch III were calculated and normalized with respect to GTP-bound Ga S .We observed a marked decrease in the nonbonded interaction energy for both GTP-bound and GDPbound Ga S [R228C] (Fig. S5), consistent with analysis of the hydrogen bond network and decreased T m value for GTPbound Ga S [R228C].Relative to Ga S [R228C], we observed an increase in the nonbonded interaction energy for both GTPbound and GDP-bound Ga S [R258A] (Fig. S5).This is also consistent with our analysis of the hydrogen bond network for this variant and the increased T m value for GDP-bound Ga S [R258A].Interestingly, we did not observe decreased interaction energies for GDP-bound Ga S or GDP-bound Ga S [Q227].Thus, while quantification of interaction energies was consistent with observed T m values for two of the variants, analysis of hydrogen bonding networks appears to be more qualitatively consistent with observed differences in protein stability.

STD-NMR spectroscopy with Ga S and Ga S variants in complexes with nucleotides
From the allosteric model connecting the nucleotidebinding pocket to changes in interactions among the switch region, we hypothesized that the Ga S disease-associated mutation should also alter interactions between nucleotides and nearby residues within the Ga S nucleotide-binding pocket, even for sites of mutations that do not directly interact with nucleotides.We utilized STD-NMR spectroscopy to investigate the local structure and dynamics of residues within the nucleotide-binding pocket for the seven disease-associated variants.STD-NMR has been employed extensively to provide information on protein-ligand interactions and is especially suited for elucidating information with ligands that exhibit moderate binding affinities (K d > 1 mM) (28)(29)(30)(31)(32). Wellisolated resonances from the protein, such as methyl groups that exhibit unique chemical shifts, are saturated in the 1 H-NMR spectrum.If a ligand interacts with the protein, polarization is transferred to the ligand from the protein via spin diffusion (29).By subtracting the saturated spectrum from a reference spectrum, a difference spectrum is obtained where signals from protons of the interacting ligand are Disease-associated mutations of G proteins predominantly observed.The intensities of the observed 1 H ligand signals depend on both the average distance and the dynamics of noncovalent interactions between the protein and ligand (29,31).
STD-NMR experiments should be sensitive to changes in the local structure and dynamics within the nucleotide-binding pocket and provide complementary information to the CD measurements and MD simulations observing changes in interaction networks.For example, while the variabletemperature CD data and MD simulations did not reveal clear-cut differences between Ga S and Ga S [Q227L], we hypothesized that local differences in the nucleotide-binding pocket may be observed by STD-NMR.Additionally, we were interested to explore whether STD-NMR data pointed to differences in potential mechanisms among the loss-of-function mutations Ga S [R258A], Ga S [R265H], and Ga S [R228C].
We recorded one-dimensional 1 H STD-NMR data with Ga S and the 7 Ga S variants in the presence of a 50-fold excess of GDP or GppNHp (see Experimental procedures for more details).GppNHp was used rather than GTPgS due to different availabilities of the two compounds at the time of the study.Both GppNHp and GTPgS are widely used in studies of G proteins as both have been shown to exhibit affinities for most G proteins that differ by less than an order of magnitude (33,34). 1 H NMR signals for both GDP and GppNHp were identified and assigned in one-dimensional reference spectra (Fig. S7) by transferring assignments reported from earlier studies (Biological Magnetic Resonance Bank accession code bmse000270) (35).From exploratory STD-NMR experiments, we identified four 1 H signals that were well resolved and free from residual protein background signals that were assigned to protons in the guanosine and sugar ribose groups, labeled as protons "a" through "d" (Fig. S7).STD-NMR experiments with GDP utilized signals from these four protons; in experiments with GppNHp, protons "c" and "d" were overlapped with background signal from the protein, and we therefore excluded them from analysis of samples containing GppNHp.Examination of the structure of Ga S (PDB ID 1AZT) showed multiple residues within $6 Å of the protons monitored in STD-NMR experiments (Fig. S7).We anticipated these residues would make the largest contributions to observed STD-NMR signals and predicted that either structural rearrangement of these residues or their fluctuation relative to the bound nucleotide would result in changes in observed STD-NMR signals.
In initial STD-NMR experiments with Ga S and GDP, we optimized the experimental signal-to-noise by varying the saturation transfer time from 0.5 to 3.0 s (Fig. S7).These data confirmed that interactions between nucleotides and Ga S were detected in STD-NMR experiments and allowed optimization of the transfer time.We observed a significant increase in the STD-NMR signal up to 2.0 s, and only minor increases in STD-NMR signals were observed with a 3.0 s transfer time.We therefore selected 2.0 s for all subsequent experiments.Control experiments recorded for samples containing either only Ga S and buffer or only nucleotide and buffer showed no STD-NMR signals, as expected (Fig. 4).This confirmed that signals observed in STD-NMR experiments were due to specific interactions between nucleotides and Ga S or Ga S variants.
The observation of STD-NMR signals for all variants provided important confirmation that nucleotides were bound to all variants.An earlier study of Ga S [R265H] reported that this variant did not bind GTP or GTP analogs (36), though this was later questioned (21).Importantly, we observed STD-NMR signals for Ga S [R265H], indicating that Ga S [R265H] does bind both GDP and GTP (Fig. 4) and is more in line with results by Hu and Shokat (21) (see Discussion).Detection of STD-NMR signals for all other Ga S variants also confirmed their abilities to interact with and bind both GDP and the GTP analog.
In STD-NMR spectra of Ga S or Ga S variants with GDP, we observed significant variation in the signal intensities of protons "a" through "d" (Figs. 4 and S8).Mini-Ga S showed the smallest signal intensities in STD-NMR spectra with respect to Ga S for protons "b" through "d," which appears consistent with its expected weaker nucleotide binding.The relatively largest signals for GDP were observed with the two variants Ga S [R258A] and Ga S [R265H] (Fig. 4).In STD-NMR spectra of Ga S or Ga S variants with GppNHp, we also observed a large variation in signal intensities, with mini-Ga S again showing the weakest signal intensities (Figs. 4 and S8), confirming that the STD-NMR data were sensitive to the degree of interaction between the ligand and protein.
We integrated the observed signal intensities from the STD-NMR spectra and used them to calculate corresponding STD amplification factors (STD-AFs) by multiplying the observed STD-NMR signal by the molar excess of ligand over protein (see Experimental procedures), following previously validated protocols (28,29).The STD-AF converts the STD signal intensity, which depends on the fraction of bound ligand, to a value that depends on the fraction of bound protein (28,29).We then calculated and compared the STD-AF values for all samples of Ga S and Ga S variants with GDP and GppNHp (Fig. 5).The STD-AF values were normalized with respect to the "a" proton of the guanine nucleobase, because in MD simulations, the immediate environment around the "a" proton showed little change (see below and Fig. S9).
The normalized STD-AF values for samples of Ga S or Ga S variants prepared with either GDP or GTP provided a "fingerprint" of the extent of interaction of the nucleotide protons with local protons from the protein.For Ga S , relative to the "a" proton, we observed an increased STD-AF for the b and c protons with GDP, and a decrease in the STD-AF for the "b" proton with GppNHp (Figs. 5 and 6).We observed striking differences in the STD-AF profiles for both GDP-bound Ga S [R258A] and Ga S [R265H], as compared with Ga S and all other Ga S variants (Figs. 5 and 6).For both Ga S [R258A] and Ga S [R265H], we observed significant increases in their STD-AF values for complexes with GDP.This appears to qualitatively correlate with the rapid dissociation rate of GDP for these variants (36).For Ga S [Q227L], we observed significantly lower STD-NMR values than Ga S .This suggests that while the global conformations and switch region interaction Disease-associated mutations of G proteins networks are similar between Ga S [Q227L] and Ga S , the local structure and dynamics of residues within the nucleotide binding are different and reflected as differences in the STD-NMR data.We observed the smallest STD-AF values with GDP for mini-Ga S .Though the STD-AF values do not necessarily directly reflect binding affinities, the lower STD-AF values observed for mini-Ga S are consistent with our expectations that this engineered protein exhibits significantly altered binding with nucleotides (24), which is also consistent a loss of interactions between bound nucleotides and residues nearby in both switch II and switch III regions removed to facilitate crystallization.
For protein complexes with GppNHp, STD-AF values calculated for "a" and "b" protons show a decrease in the value of the "b" proton relative to "a" for Ga S (Figs. 5 and 6), establishing a "fingerprint" of interaction for the native protein.
Compared with the STD-AF profile of Ga S , the STD-AF values for five of the Ga S showed similar profiles, with two exceptions.Both Ga S [R258A] and Ga S [R228C] showed larger values observed for the "b" proton for their complexes with GppNHp (Figs. 5 and 6).While we did observe an STD-NMR signal for GppNHp in the presence of mini-Ga S , the value for the "b" proton was the lowest observed within the study, consistent with the expectation of altered binding of this engineered protein to the GTP analog (Figs. 5 and 6).
Observations from the STD-NMR data reflect changes in the local structure and dynamics of amino acids within the nucleotide-binding pocket of Ga S and Ga S variants; however, different configurations of the nucleotide binding region may result in similar STD-AF values.To obtain a more detailed view, we analyzed our MD simulations to investigate potential differences in the dynamics and configuration of the binding pocket, focusing on GDP-bound and GTP-bound Ga S and the three variants explored in the earlier structural investigations Reference" is a one-dimensional 1 H NMR spectrum of GDP, "protein only" is an STD-NMR control experiment with a sample containing Ga S and buffer but no nucleotide, and "nucleotide only" is a STD-NMR control experiment with a sample containing nucleotide and buffer but no protein.B, the chemical structure of GppNHp and one-dimensional STD-NMR spectra of complexes with GppNHp.Other presentation details are the same as in A. Views of the presented NMR data were expanded from the full spectra shown in Fig. S7.STD, saturation-transfer difference.

Disease-associated mutations of G proteins
underlying different thermal melting behavior, Ga S [R228C], Ga S [Q227L], and Ga S [R258A].For each protein and for each nucleotide complex, we monitored each of the four protons on the nucleotide, "a" through "d," observed in the STD-NMR experiments over the course of the simulations.In each simulation, we identified the closest protons in the protein within a radius of $6 Å of the nucleotide and then calculated the fraction of the simulation time that each residue spends in proximity to one or more of the monitored protons (Figs. 7, S9 and S10).
We observed both similarities and differences in the local environment among Ga S and the 3 Ga S variants.The "a" proton on the guanine ring of the nucleotide is in persistent proximity with G52 for both GDP-bound and GTP-bound Ga S and for all variants, with an overall mean occupancy across all proteins of 97.4 ± 1.4% (Figs. 7 and S9).The associated SEM was used as an indicator of the variability, which was relatively low in this case.As discussed above, this mostly uniform local environment around proton "a" rationalized our normalization of the STD-AF values with respect to its signal intensity (Figs. 5 and 6).
The "b" proton is mainly interacting with Asp173 (mean occupancy of 67 ± 14%) and Lys293 (32.2 ± 14.1%, Fig. S9).Most of this variability is due to the significant difference between the GDP and GTP-bound forms in both Ga S [Q227L] and Ga S [R258A] mutants, where the GTP-bound form shows an even contact between the "b" proton and Asp173 and Lys293.Interestingly, the STD-AF normalized values of these two mutants show a gap for the "b" proton interactions between their corresponding GDP and GTP complexes.For the remaining protons, we were unable to experimentally measure their GTP-bound normalized STD-AF values.Nevertheless, we can examine results from MD simulations to make predictions on their behavior.Proton "c" exhibits the highest variability in the identity of the closest residue among the different Ga S forms with shared interactions with Arg201, Asp173, and Glu50; however, for Ga S [Q227L], proton "c" almost exclusively interacts with Glu50 when bound to GDP or with Arg201 when bound to GTP.Analysis of the "d" protons in GDP was more complex since this signal is due to two equivalent protons.To address this, we analyzed the arithmetic average of the independent measurements for each proton, and the results show a constant and even distribution of the closest interacting residues, with Glu50 and Arg201 making an even interaction with the "d" protons in all simulations.Moreover, in all four protein systems (Ga S or variants Ga S [Q227L], Ga S [R258A], and Ga S [R228C]) the GTP-bound form displays a distinct interaction of the "d" protons with Figure 6.STD-NMR amplification factors mapped onto the chemical structures of GDP and GppNHp for complexes with Ga S and Ga S variants.For each proton observed in the STD-NMR spectra, the relative STD-AFs are colored as a percentage of the STD-AF normalized with respect to the "a" proton (same as in Figs. 3 and 4) indicated with the black circle.The STD-AFs are colored according to their relative intensities compared to those for the "a" proton: relatively weaker interactions representing <50% (gray), weak-moderate interactions representing 51% to 99% (yellow), moderate interactions representing 101% to 150% (green), and stronger interactions representing >150% (red).The white circles indicate that no analysis was performed for those protons due to overlap of nucleotide and background residual protein signals.AF, amplification factor; STD, saturation-transfer difference.
Gly52 (average occupancy 13.0 ± 3.8%, versus 0.6 ± 0.8% in the corresponding GDP-bound forms, Fig. S9), suggesting a different microenvironment of these protons as a function of the nucleotide bound.

Discussion
Based on the reported crystal structure of Ga S [R201C] and biochemical data of Ga S [R228C], mechanisms for Ga S gain-offunction and loss-of-function mutations were proposed that related changes in hydrogen bond networks within the switch regions to changes in protein stability and function (21).This framework was built upon earlier biophysical and biochemical experiments with Ga i1 that related changes in Ga i1 thermal stability to changes in the hydrogen bond network among switch region residues, which were also closely related to global structural rearrangements and changes in Ga i1 function (18,19,26).We leveraged an integrative biophysical approach to investigate the extent to which this framework could be expanded to provide mechanistic insights into additional disease-associated Ga S variants.
With this biophysical toolbox, we examined four loss-offunction Ga S variants: Ga S [G226A], Ga S [R228C], Ga S [R258H], and Ga S [R265H].Based on biochemical data, it had been proposed that the R228C mutation resulted in loss of function due to destabilization of its GTP-bound conformation (21).We observed a decreased T m for GTP-bound Ga S [R228C] with respect to GTP-bound Ga S , indicating decreased stability (Fig. 2), and we correlated this observation with the loss of hydrogen bonds in the switch region network for this variant (Fig. 3).We also observed the largest RMSFs within the switch region for GTP-bound Ga S [R228C] (Fig. S4).Taken together, these data are consistent with the explanation of a destabilized GTP-bound state for this variant.Our observations of decreased stability for GTP-bound Ga S [G226A], closely similar to Ga S [R228C], suggest that the same mechanism likely explains the impact of the G226A variant as well.
For the loss-of-function variants Ga S [R258A] and Ga S [R265H], we observed the stability for these GTP-bound variants to be similar to that of Ga S , in contrast to Ga S [G226A] and Ga S [R228C] (Fig. 2).For Ga S [R258A] and Ga S [R265H], we also observed increased stability for their GDP-bound states, which was found to be consistent with increased hydrogen bond networks in MD simulations (Figs. 2  and 3).While our observations for Ga S [G226A] and Ga S [R228C] confirmed the correlation between increased stability and increased interactions within the switch regions, this was not consistent with the behavior of other loss-offunction variants Ga S [R228C] and Ga S [G226A], indicating that an alternative mechanism must account for their loss of function.
Interestingly, for nucleotide free Ga S [R258A] and Ga S [R265H], we observed T m values that were similar to the higher T m value observed for GTP-bound Ga S (Fig. S3).The increased stability of nucleotide-free Ga S [R258A] and Ga S [R265H] indicates that these two variants adopt stable, relatively rigid states in the absence of any nucleotide.This appears to be a unique property of these variants and is not shared with nucleotide free Ga S or Ga S [R201C].Based on these data, we propose an alternative mechanism to explain the loss of function in Ga S [R258A] and Ga S [R265H].Specifically, we propose that Ga S [R258A] and Ga S [R265H] are "locked" in a conformation that is suboptimal for binding GDP or GTP, impairing the protein's ability to undergo conformational changes that support stronger interactions with nucleotides, resulting in rapid dissociation of GDP and rapid hydrolysis of GTP.This alternative mechanism is supported by earlier studies that reported Ga S [R265H] did not bind GTP (36).Leyme et al. based their conclusion that Ga S [R265H] did not bind GTP based on the lack of increased inherent tryptophan fluorescence upon the addition of GTP (36).However, more recent studies by Hu and Shokat, using different approaches that did not rely on inherent tryptophan fluorescence, concluded that Ga S [R265H] does bind GTP (21).From the STD-NMR data in the current work, we observed that Ga S [R265H] and Ga S [R258A] bind both GDP and GTP analogs (Figs.4-6).This discrepancy with Leyme et al.'s conclusion can be explained by the idea that Ga S [R265H] does not undergo conformational changes upon nucleotide exchange that would affect the environment around tryptophan residues, particularly W234, which is adjacent to the switch regions, leading to no change in tryptophan fluorescence intensity.One implication from this idea is that the rigidity of Ga S [R265H] and Ga S [R258A], even in the nucleotide-free state, suggests these variants may be less suitable as targets for potential small The color palette for the residues interacting with these protons denotes increasing shades of pink which are proportional to the frequency of the residue as being the closest to the given proton, with darker shades of pink indicating higher frequencies as a function of the mutations explored.STD, saturation-transfer difference.

Disease-associated mutations of G proteins
molecule therapeutics that require access to flexible binding pockets.
For the gain-of-function mutations Ga S [R201C], Ga S [Q227L], and Ga S [A366S], based on the earlier crystal structure and mechanistic study of GDP-bound Ga S [R201C] (21), we wanted to first test to what extent their relative stabilities correlated with changes in hydrogen bonding networks, and, second, to what extent these observations could be correlated with their changes in functional activity.For Ga S [R201C], we observed an increased stability and increased hydrogen bond network (Figs. 2 and 3) for GDP-bound Ga S [R201C].This is consistent with the earlier claim of increased stability for this variant based on the crystal structure (21).For Ga S [A366S], we observed a lower T m for the GTP-bound conformation, indicating a loss of stability and likely loss of hydrogen bonding network for GTP-bound Ga S [A366S], destabilizing the GTP-bound conformation of this variant (Figs. 2 and 3).Earlier measurements of Ga S [A366S] nucleotide binding kinetics showed a decreased ability of this variant to bind GDP (37), suggesting a plausible contribution toward its increase in function.Our analysis indicates that another contributing factor toward increased function is the destabilization of the GTP-bound conformation.Though this observation differs from the increased stability observed for GDP-bound Ga S [R201C], it is line with earlier studies of A326S in Ga i that showed this mutation altered the conformation of the nucleotide-binding pocket through steric crowding (38), which is consistent with our observations of a destabilized GTP-bound conformation.
Surprisingly, Ga S [Q227L] is the only variant that exhibited stability for both GDP and GTPgS complexes highly similar to native Ga S (Fig. 2).This was the only variant that did not exhibit a change in stability for either the GDP-bound or GTPbound states.Consistent with this result, we also observed similar hydrogen bonding networks between Ga S [Q227L] and Ga S (Fig. 3), indicating that this particular mutation does not inhibit this variant's ability to adopt different global conformations for complexes with different nucleotides.Several unique features were observed in MD simulations for Ga S [Q227L] (Figs. 7 and S9), providing potential clues to its gain in function.For Ga S [Q227L], proton "c," located on the pentose sugar ribose of GDP, exhibited significant differences in its local environment as compared with Ga S , Ga S [R228C], and Ga S [R258A].Within the local environment of proton "c," for Ga S [Q227L], we observed a significant increase in the percent contact with residue Glu50 (96.3%) and loss of contact with residues Asp173 and Arg201.While these differences are more subtle than for other variants, they may partially contribute toward functional differences and reflect differences in the dynamics of residues within the nucleotide-binding pocket between Ga S [Q227L] and Ga S .
In summary, from our combined biophysical and computational data, we observed a clear correlation between protein stability and changes in hydrogen bonding within the switch region for a series of seven disease-associated Ga S variants containing single amino acid mutations.Six of the seven variants could be identified and categorized based on observations of their stabilities, suggesting that experiments measuring the stability of additional Ga variants may be sufficient for providing a plausible mechanism of loss or gain of activity.Comparison of the relative stabilities of these variants for complexes with nucleotides, the hydrogen bonding network within the switch regions, and structural dynamics within the nucleotide-binding pocket provided evidence supporting five mechanisms that give rise to changes in function.Four of these mechanisms are related by the observation that changes in stability for either GDP-bound or GTP-bound variants correlate with changes in hydrogen bond networks and can rationalize protein function.In the first two related mechanisms, Ga S [R201C] results in a gain of function by stabilization of the GDP-bound conformation while Ga S [R228C] and Ga S [G226A] result in loss of function by destabilization of their GTP-bound conformations.In a third mechanism, Ga S [R265H] and Ga S [R258A] result in a loss of function by being locked into conformations that are not optimal for binding GDP and that also result in rapid GTP hydrolysis.In a fourth mechanism, Ga S [A366S] results in destabilization of the GTP-bound conformation, possibly due to steric crowding around the nucleotide-binding pocket.Evidence for these four mechanisms could be observed from the combined CD and NMR data and MD simulations relating changes in variant stability with changes in hydrogen bond networks and function.The variant Ga S [Q227L] was the only one that could not be readily distinguished from the stability or hydrogen bonding interactions observed for Ga S , suggesting a more subtle mechanism underlies its activity.We anticipate that by making similar measurements with additional Ga variants, future research may be able to categorize mechanisms of loss or gain of function for many other variants based on similar biophysical data.

Generation of Ga S and variant Ga S plasmids
An initial Ga S construct was gifted from the laboratory of Prof. Kevan Shokat (UCSF) (21), contained in a pET15b vector and derived from the short human isoform of Ga S (PubMed accession number NP_536351) with the first seven residues deleted.We modified this plasmid by replacing a DRICE cleavage site with a tobacco etch virus nuclear-inclusion-a endopeptidase (TEV protease) cleavage site.This plasmid was used to express Ga S for our experiments.All Ga S variants were obtained from this plasmid by site-directed mutagenesis using the QuikChange II kit (Agilent) and primers listed in Table S2.The plasmid encoding for mini-Ga S was obtained from GenScript and designed to be identical to the sequence of the mini-Ga S "393" sequence (24).Disease-associated mutations of G proteins 34 mg/ml chloramphenicol and 100 mg/ml carbenicillin.Cells were grown to an absorbance (A 600 nm) between 0.6 and 0.8 and protein expression was then induced with 50 mM IPTG at 25 C for 18 to 20 h.Cells were sedimented by centrifugation at 4000g at 4 C for 20 min then were lysed using a cell disruptor at 24,000 psi and 4 C in buffer containing 50 ml of 25 mM Tris pH 8.0, 150 mM NaCl, 1 mM MgCl 2 , 50 mM GDP with EDTA-free in-house protease inhibitor cocktail (500 mM

Expression and purification of mini-GaS
We expressed and purified mini-Ga S using a strategy adapted from earlier studies (24,27).BL21-CodonPlus (DE3)-RIL were transformed with the expression vectors and grown at 30 C in Terrific broth medium with 0.2% glucose, 34 mg/ml chloramphenicol, 100 mg/ml carbenicillin, and 5 mM MgSO 4 .Protein expression was induced with 50 mM IPTG at 25 C and left to grow overnight (16-20 h).Cells were lysed in 40 mM Hepes pH 7.5, 100 mM NaCl, 10 mM imidazole, 10% vol/vol glycerol, 5 mM MgCl 2 , 50 mM GDP with EDTA free in-house protease inhibitor cocktail (500 mM 4-(2-aminoethyl)-benzenesulfonylfluoride hydrochloride, 1 mM E-64, 1 mM leupeptin, 150 nM aprotinin), and using a cell disruptor operating at 24 kPsi and 4 C.The lysate was clarified by centrifugation at 25,000 rpm and 10 C for 30 min and loaded on a His-Trap HP NiNTA 5 ml column (Cytiva) via an Akta Start FPLC.Resin was washed with lysate buffer containing a final concentration of 40 mM imidazole and 500 mM NaCl, and protein was eluted from the resin with lysate buffer containing 400 mM imidazole.The sample was exchanged into buffer via an Akta Start FPLC equipped with a HiPrep 26/10 desalting column (Cytiva) that was equilibrated in 20 mM Hepes pH 7.5, 100 mM NaCl, 10% vol/vol glycerol, 1 mM MgCl 2 , and 10 mM GDP. Fractions containing mini-Ga S were collected and incubated for 16 to 20 h at 4 C with thermostabilized TEV protease at 1:100 M ratio (mini-Ga S :TEV).Cleaved mini-Ga S was purified via reverse immobilized mobile affinity chromatography by incubating the sample with nickel resin for 30 min and collecting the flow through using a gravity flow column.Mini-Ga S was further purified via an Akta Start FPLC equipped with a HiLoad Superdex 75 pg column (Cytiva) equilibrated in storage buffer (10 mM Hepes pH 7.5, 100 mM NaCl, 10% vol/vol glycerol, 1 mM MgCl 2 , 1 mM GDP, and 0.1 mM tris(2-carboxyethyl)phosphine).Purified samples were rapidly frozen in liquid nitrogen and stored at −80 C until needed for experiments.

Tryptophan fluorescence assays
Ga S and Ga S variants were exchanged into buffer containing 20 mM Hepes pH 8.0, 150 mM NaCl, and 5 mM GDP using a PD MiniTrap G25 column (Cytiva), and the sample volumes were adjusted to achieve a final concentration of 10 mM protein.Tryptophan fluorescence spectra were recorded with a Cary Eclipse spectrometer operating at 25 C.The excitation wavelength was 280 nm and emission was observed at 340 nm with 5 nm excitation and emission bandwidths.GTPgS was added to 200 ml protein for a final concentration of 50 mM in a quartz cuvette.The sample was mixed vigorously by repeated pipetting and data collection immediately started.Data points were recorded every 12 s for 200 min total acquisition time.Data were analyzed in GraphPad prism (https://www.graphpad.com)using the one-phase association equation shown below: where Y o is the fluorescence intensity when x (time) equals zero, "Plateau" is the fluorescence intensity when the system has been fully equilibrated at infinite time expressed in the Disease-associated mutations of G proteins same units as Y , and k is the rate constant, expressed as the reciprocal of the x-axis time units.

Protein thermal unfolding monitored by CD
Single-wavelength variable temperature CD data were recorded on a Chirascan qCD spectrometer.Ga S or Ga S variants were concentrated to 10 mM in buffer containing 10 mM Hepes pH 7.5, 100 mM NaCl, 1 mM MgCl 2 , and 50 mM GDP or 50 mM GTPgS.First, full-wavelength CD spectra were recorded at 25 C to confirm samples were properly folded.Then, protein melting was monitored by recording the ellipticity at 220 nm from 25 C to 89 C with a linear heating rate of 1 C/min.Each T m was determined by fitting the data in GraphPad prism to a Boltzmann sigmoidal equation shown below: where Y is the fraction of folded protein as a function of temperature, V 50 is the temperature at which half of the protein is unfolded, X represents the temperature, and "slope" is the steepness of the curve."Top" and 'Bottom' are normalization factors set to 1.0 and 0, respectively.Error bars were determined by calculating the standard deviation from three independent measurements.

NMR experiments
Samples of Ga S or Ga S variants were exchanged into NMR buffer using a PD MiniTrap G25 column (Cytiva) equilibrated with 25 mM HEPES pH 7.0, 75 mM NaCl, 5 mM MgCl 2 , 100 mM DSS, 10% v/v 2 H 2 O, and either 2 mM GDP or 2 mM GppNHp.Exploratory experiments employing a range of protein-to-ligand concentrations were tested to determine optimal STD-NMR signal-to-noise.Based on these data we prepared samples containing 40 mM protein and 2 mM ligands (Fig. S6).NMR experiments were recorded at 25 C and 800 MHz 1 H Larmor frequency using a Bruker Avance III spectrometer equipped with a 5 mm TXI cryoprobe.Onedimensional-1 H saturation transfer difference (STD) experiments were recorded using water suppression by gradient-based excitation sculpting (Bruker pulse program stddiffesgp), with 64 scans per experiment and a 3 s delay between scans.Onresonance spectra were acquired by irradiating at 0.79 ppm to target Ga S methyl protons, and off-resonance spectra were acquired by irradiating at −40.0 ppm.Saturation times were tested from 0.5 to 3 s, and an optimal saturation time of 2 s was determined (Fig. S8).Data were acquired and analyzed in Topspin 4.1.1.Prior to Fourier transformation, the data set was multiplied by an exponential window function applying 3 Hz line broadening and zero-filled to 65,536 points; baseline correction was applied to all frequency domain data sets.To reduce background signals, a reference STD-NMR spectrum was recorded with a sample containing the same concentration of protein but no nucleotide.This spectrum was processed identically and subtracted from all STD-NMR spectra.
AFs, A STD , were determined from the integrals of ligand signals in the STD difference (I 0 − I sat ) and reference (I 0 ) spectra according to the following equation from reference (29) where A STD is the AF of each STD signal, I O is the integral of the STD signal in the reference spectrum, I sat is the integral of the STD signal in the difference spectrum and [L] and [P] are the concentrations of ligand and protein, respectively.Signal-to-noise ratios in the STD difference spectra, SNR diff , were determined for each ligand signal using the "sinocal" function on topspin 4.1.1.Error associated with the AFs, s ,was evaluated using the following equation: where A STD is the AF of each STD signal and SNR diff is the signal-to-noise ratio of the STD signal.For presentation of the normalized A STD signals, each nucleotide signal "a" to "d" and associated errors were normalized by the A STD of the "a" signal.

MD simulations
The crystal structure of the Ga S protein in complex with GTPgS was retrieved from the Protein Data Bank (PDB ID 1AZT) (40) and prepared for MD simulations as follows.Chain B from the dimer in the asymmetric unit was retained, and the missing segment between a1 and aA (residues 70-86) were modeled with Prime (Schrӧdinger Release 2011-2: Prime; Schrӧdinger, LLC: New York, 2011).Crystallographic water molecules beyond 5 Å of GTPgS as well as all cocrystallized phosphate ions were removed, and the protonation states of ionizable residues were determined with PROPKA at pH 7.0.The GTPgS substrate was modified to either GTP or GDP, depending on the intended system, with Maestro v 9.4 (Schrӧdinger, LLC).Amino acid replacements to generate the different Ga S variants were produced with the Schrödinger's PyMOL mutagenesis tool (The PyMOL Molecular Graphics System, Version 2.0.Schrӧdinger, LLC.).
MD simulations were executed with the Q software package (41), using the OPLS-AA/M force field for proteins and the TIP3 water model (42), available from https://github.com/qusers/q6.The parameters to describe GTP and GDP were taken from same force field.All models were partially solvated in a water sphere of diameter 50 Å with its center located at the nucleotide's center of geometry.Protein atoms outside the sphere were tightly constrained to their initial coordinates with a force constant of 200 kcal mol −1 Å −1 and excluded from nonbonded interactions, while water molecules at the sphere boundary were subjected to radial and polarization restrains following the SCAAS model (43).All atoms inside the simulation sphere were allowed to move freely, where the local Disease-associated mutations of G proteins reaction field multipole expansion was used for long-range electrostatic interactions beyond a cut-off of 10 Å (44).Each system was subject to 100 independent replica MD simulations, each of them starting with an equilibration phase consisting of 51 ps of gradual heating to 298 K and concurrent release of (initial 10.0 kcal mol −1 Å −2 ) harmonic restraints on solute heavy atoms.Unrestrained MD simulations followed for 5 ns at a temperature of 298 K and employed a 1 fs time step, giving a total of 0.5 ms of unrestrained simulation per system, accumulating 2 ms MD simulations in this study.

Figure 1 .
Figure 1.Crystal structure of Ga S with annotated disease-associated mutations.A, the GTP-bound Ga S crystal structure (PDB: 1AZT) (40).GTP is shown in stick representation and disease-associated mutations shown to increase or decrease adenylyl cyclase activity is shown in cyan or red stick representation, respectively.G226 is shown as a cyan circle.Magnesium is shown as a green sphere.B, expanded view of the nucleotide-binding pocket and annotated mutations of interest.Same presentation details as A.

Figure 2 .
Figure2.Stability of GaS and GaS variants determined by variable-temperature CD spectroscopy.A, representative CD spectrum of GaS bound to GDP (blue curve) or GTPgS (orange curve).B, thermal unfolding of GaS bound to GDP or GTPgS monitored by variable temperature single wavelength CD spectroscopy.Same color scheme as A. C, histogram of Tm values determined from fitting the variable temperature CD data of GaS, mini-GaS, and GaS variants obtained in complex with GDP or GTPgS.Variant names are organized according to those with predicted increased stability for the complex with GDP (cyan) or predicted decreased stability for the complex with GTP (red), relative to GaS.Dashed horizontal lines are located at the mean Tm values for GaS in complexes with GDP or GTPgS.Error bars were determined from the SD of triplicate measurements.

Figure 3 .
Figure 3. Analysis of hydrogen bond interactions between switch regions extracted from MD simulations.A, hydrogen-bond interaction profiles for Ga S (labeled "WT"), Ga S [Q227L], Ga S [R228C], and Ga S [R258A], focused on interactions between switch I, switch II, and switch III regions.Colored boxes indicate the existence of a persistent hydrogen bond, determined with different thresholds of occurrence along the given MD simulation, proportional to the color intensity for GDP-bound (blue) or GTP-bound (orange) proteins.B, schematic representation of all hydrogen-bond interactions depicted in the Hbond matrices.C, the structure of GTP-bound Ga S with the involved residues from the analysis annotated.Residues from switch I are colored yellow, switch II colored blue, and switch III colored red.MD, molecular dynamics.

Figure 4 .
Figure 4. One-dimensional 1 H STD-NMR spectra of Ga S and Ga S variants in complexes with GDP or GppNHp.A, the chemical structure of GDP and one-dimensional STD-NMR spectra of complexes with GDP.Protons observed in the STD-NMR spectra are annotated "a" to "d." Below the chemical structure are expanded regions from one-dimensional STD-NMR spectra shown in Fig. S8 containing the annotated signals."Reference" is a one-dimensional 1 H NMR spectrum of GDP, "protein only" is an STD-NMR control experiment with a sample containing Ga S and buffer but no nucleotide, and "nucleotide only" is a STD-NMR control experiment with a sample containing nucleotide and buffer but no protein.B, the chemical structure of GppNHp and one-dimensional STD-NMR spectra of complexes with GppNHp.Other presentation details are the same as in A. Views of the presented NMR data were expanded from the full spectra shown in Fig.S7.STD, saturation-transfer difference.

Figure 5 .
Figure 5. Normalized STD amplification factors measured for GDP and GppNHp in complex with Ga S or Ga S variants.Normalized STD amplification factors (STD-AF) determined for GDP (blue) and GppNHp (orange) in complex with Ga S or Ga S variants.The specific protons used to quantify the STD-AF values are annotated on the chemical structure of GDP.STD amplification factors were normalized to the "a" proton.Error bars were determined from the signal-to-noise ratios in the corresponding NMR spectrum.STD, saturation-transfer difference; AF, amplification factor.

Figure 7 .
Figure 7. Variation in the structure and dynamics of the nucleotidebinding pocket among Ga S and Ga S variants.Nucleotide-binding pocket of Ga S bound to GTP s(shown in orange stick representation, PDB:1AZT) with protons seen in STD-NMR interactions annotated "a" through "d."The color palette for the residues interacting with these protons denotes increasing shades of pink which are proportional to the frequency of the residue as being the closest to the given proton, with darker shades of pink indicating higher frequencies as a function of the mutations explored.STD, saturation-transfer difference.

Table 1
Properties and predicted impacts of Gas disease-associated mutations